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I. INTRODUCTION 



There appear to be deep connections between phenomena at the microscopic quantum level and those at the 
macroscopic level described by general relativity. These connections were first elucidated by studies of quantum field 
theory in the classical curved spacetime of general relativity, particularly those involving particle creation by strong 
gravitational fields, such as exist in the early universe and near black holes Jl],||]. The framework of quantum field 
theory in the classical curved spacetime of general relativity appears to be quite robust, and is the one we employ 
in our present study. We are particularly concerned with certain non-perturbative gravitational contributions to the 
vacuum energy of quantized fields j| . 

We show that these contributions may account for the recent observations |^|| of type-la supernovae (SNe-Ia) 
that seem to imply that there is an acceleration of the recent expansion of the universe. In our model, the key 
new ingredient needed to account for the observations is the existence of a particle having a very small mass of 
about 10~ 33 eV. This could be a scalar particle or one of higher spin, such as a graviton. It is well-known that in 
a Robertson- Walker universe, the linearized Einstein equations obeyed by the graviton field take the same form (for 
each polarization in the Lifshitz gauge) as the equation of a minimally coupled massless scalar field . Inclusion of 
a mass of about 10 -33 eV in these equations would give effects similar to the ones studied here, and would evidently 
not conflict with other observations. We consider a scalar field for simplicity. 

The effects we investigate stem from the discovery |3|,|| that a covariant infinite series of terms in the propagators 
of quantum fields in curved spacetime can be summed in closed form, to all orders in the curvature. The summed 
infinite series of terms are all those that involve at least one factor of the scalar curvature, R, (together with any 
number of factors of the Riemann tensor and its covariant derivatives) in the Schwinger-DeWitt proper-time series 
for the propagator. We will refer to this partially-summed form of the propagator as the R-summed form of the 
propagator. Although only the first several terms in the proper-time series have been calculated (because of their 
complexity), the R-summed form can be derived by means of mathematical induction A further reason to expect 
the R-summed form of the propagator to contain information of physical significance, is that the leading term in 
the R-summed form of the propagator follows directly from the Feynman path-integral expression for the propagator 
by means of a Gaussian integration about the dominant path || . Our results flow from the effective action that is 
obtained from the R-summed form of the propagator. 

We leave for a later paper, the contribution of these non-perturbative gravitational effects to early inflation at 
the grand-unified scale and associated particle creation. In this paper, we focus on possible consequences of these 
non-perturbative terms in the effective action that may be observed today. As noted earlier, we find that these non- 
perturbative effects in the presence of an ultralight particle having a mass of about 10 -33 eV, give a good fit to the 
observed SNe-Ia data points. The effective action of our model has a single free parameter determined by the mass of 
the particle and its curvature coupling constant. In addition, there are two more parameters which characterize the 
solutions to the effective gravitational field equations. These can be taken to be the present Hubble constant and the 
present matter density. The Hubble constant is determined by low redshift measurements as usual. The mass of the 
particle, and the present matter density are given by a fit to the SNe-Ia data. The age of the universe, and the cosmic 
time, or red-shift, at which the non-perturbative contributions of the ultralight particle first become significant, follow 
from these parameters. 

Our model involves a renormalizable free field theory in curved spacetime, and apart from the supposed existence of 
an ultralight mass, is based on previously discovered non-perturbative terms. Furthermore, it requires no cosmological 
constant in the usual sense of the term; and no fine-tuning of the value of the current cosmic time is necessary to 
explain the current observed value of the matter density. 

The outline of this paper is as follows. In Sec. II, we calculate, by zeta-function regularization, the renormalized 
effective action that follows from the R-summed propagator. This confirms the result found with dimensional regular- 
ization by Parker and Toms H , and generalizes it to include the case when there is an imaginary part of the effective 
action. Variation of this effective action gives the Einstein gravitational field equations, including vacuum contribu- 
tions of quantized fields. In addition, an imaginary part of the effective action implies a rate of particle creation in 
the same way as it does in quantum electrodynamics jlCfl . The effective action has infrared- type divergences which 
are treated separately in Appendix A. 

The effective action at low curvatures gives rise to induced gravitational coupling constants, as is well-known. These 
constants depend on a renormalization scale parameter. In Section III, we consider the dependence of the induced 
gravitational constants on the renormalization scale, and identify the renormalized constants with their known values 
at low curvature. This identification is used to fix all the coupling constants in the effective action, in terms of their 
low curvature values, except the field mass m and its curvature coupling £. In this context, the conclusions of Ref. 

are explored in greater detail and generalized here. Variations of various terms in the effective action are listed in 
Appendix B. 
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In Section IV, we clarify the meaning of exact and perturbative solutions to the semiclassical Einstein equations. In 
particular, we emphasize that a perturbative treatment in h must be done in such a way that only genuine quantum 
corrections are treated in a perturbative manner, while factors of Ti which occur in the "classical" Klein-Gordon action 
are not treated perturbatively. 

In Section V, we display a set of constant curvature (de Sitter) solutions^, which exist even in the absence of an 
explicit cosmological constant term in the effective action. We show that, for £ — 1/6 > 0, such solutions typically have 
Planckian curvatures and are therefore not of much physical interest. On the other hand, for £ — 1/6 < 0, we show that 
there exist de Sitter solutions with scalar curvature approximately equal to m 2 / (1/6 — £), and that such solutions exist 
for a wide range of values of to and £. As argued in later sections, these solutions could play a vital role at the present 
time, if there exists an ultralight particle with £ — 1/6 < 0. In Section VI, we carry out perturbative calculations 
of the scalar curvature about a classical Robertson- Walker universe, containing classical matter and radiation. We 
find that quantum corrections to the classical scalar curvature of the universe can become significant at an epoch 
determined by the mass scale to = mj \Jl/Q — £, for £— 1/6 < 0. Furthermore, we find that the quantum corrections 
tend to drive the scalar curvature to a constant value. These considerations lead to a scenario, presented in Section 
VII, in which quantum corrections to a matter-dominated universe cause a transition to a de Sitter solution of the 
type discussed in Section V. The various cosmological parameters in such a model universe are expressed in terms 
of three basic parameters to be determined by observation: the mass scale to, the present ratio of matter density to 
critical density Oo, and the present Hubble constant Hq. 

In Section VIII, we develop this scenario further, using it to obtain magnitude-redshift curves, and comparing these 
curves to recent data from high-redshift Type la supernovae ||. The mass scale to is determined to be roughly equal 
to 10 -33 eV, while f^o is found to be less than 0.4. In Section IX, we derive from these parameters the age of the 
universe in our model, which is found to be greater than about 13 Gyr. Finally, in Section IX, we conclude with a 
discussion of our results and comment on future work. 

The main results of this paper are summarized in Eqs. ( pSj ) and (^3|) (the effective action of the theory in terms 
of measurable gravitational coupling constants) , Eq. (|5^) (the relations hip b etween scalar curvature and mass during 
the late de Sitter phase), Eq. (|8q ) (the scale factor in our model), Eq. (109) and the equations prior to it in Section 
VIII (the luminosity-distance- redshift relation), the discussion of the age of the universe in Section IX, and Figs. 3 
and 4. 

Throughout this paper, we use the metric signature convention ( — h ++), and the convention for the Riemann 
curvature tensor R^f = V WiV - T\ p ^ + T a w f av - Y a vp Y° ail . 



II. QUANTUM CORRECTIONS TO EFFECTIVE ACTION 

Here we derive the regularized effective action based on the propagator for a scalar field in curved space. 

Parker and Toms |^| define the effective action via dimensional regularization of the Feynman Green function and 
by integrating the Green function with respect to the square of the mass. Here we will define it by zeta-function 
regularization. 

To this end, consider a D-dimensional scalar field theory in curved spacetime, with classical action 

S = ~j d D x^x)H{x)cS>{x), (1) 

where H(x) = —g^{x)V^.V^. + to 2 + £R(x), and we have ignored boundary contributions. The one-loop effective 
action, W 1 , is defined by the functional integral 

e lWl ee Z = J D(j)e tS , (2) 

and is related to the transition amplitude for evolution from the "in" vacuum state at early times to the "out" vacuum 
state at late times |io[| : 

e lW± = (0„ut | in ). (3) 



Me Sitter solutions to the effective gravitational field equations have been found by other authors in gravity theories with 
higher derivative terms [^,^6). These solutions are valid at high curvatures. In this paper, we also find low-curvature de Sitter 
solutions which could play an important role in the late evolution of the universe. 
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For free fields, as considered here, the one-loop effective action, W 1 , gives the full quantum contributions of the matter 
field to the effective action. By formally performing the functional integral in Eq. (0), we get (see, for example 

W 1 = ^lndet(H/fi 2 ), (4) 

where we have introduced an arbitrary real quantity /j, with the dimensions of mass in order to render Z dimensionless. 
In terms of the zeta function, formally defined as 

((v)=TrH-», (5) 

the one-loop effective action can be expressed in the form 

W i = -l[C'(0) + ln( M 2 )C(0)]. (6) 
The proper-time heat kernel, K, is defined to satisfy the equation 

i-^-H(x))K(x,x',is) = 0, (7) 

with initial condition 

K(x.x',0) = (-g)-h^ D \x-x'). (8) 

Because Eq. (0) is a Schrodinger-like equation, the heat kernel may be formally represented as 

K{x,x',is) = (x | e~ lsIi | x'), (9) 

with the inner product (x | x') defined by Eq. (JsJ) - 

The heat kernel gives a representation for the Feynman Green's function of the theory: 



/>oo 

G(x,x') = (x I H~ l \x')= / ids{x | e -«("-« ) | x '), 

Jo 



(10) 



where we have added a small imaginary part to H for reasons of convergence. 

Also, using the Mellin transform to express the zeta function in terms of the heat kernel, one obtains 



= r(zy)- 1 / d D x^~g \ ids(i S y- x (x I e - ls(H - K) \ x). (11) 



We now introduce the i?-summed form of the heat kernel J3|,|| 

(x | e- 4s(ff " le) | x') = i(Ams)- D l 2 ^ (x, x')e l *-^ M2 - ie) F(x, x'\ is), (12) 

where A(x,x') is the Van Vleck-Morette determinant and a(x,x') is one-half the square of the geodesic distance 
between x and x' . The quantity M 2 is, in general, a function of x and x' with the property that it reduces in the 
coincidence limit to 

M 2 (x,x) = m 2 + (£-l/6)R(x) = m 2 +£R(x), (13) 

where £ = £ — 1/6. The form (|l2] ) of the heat kernel sums all terms explicitly involving at least one factor of the 
scalar curvature R, in the coincidence limit x' — > x. Away from the coincidence limit, M 2 can be taken to be the 
linear combination || 

M 2 (x, x') = m 2 + £(aR(x) + (1 - a)R(x')) (14) 

with an arbitrary choice of a, and the function F will correspondingly depend on the value of a. F may be expanded 
in an asymptotic series in powers of s, namely, 

oo 

F(x,x';is)nJ2( is yM x > x, )> ( 15 ) 
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In the coincidence limit x' — > x, F does not depend on a. We then have fo(x, x) = 1, fi(x, x) = 0, and 

J 2 (x,x) = U±-A DR+ -L(R a0lS R<*W ~ RapK*). (16) 

The fj for all j contain no factor of R (with no derivatives acting on it). 
Thus Eq. ( pi] ) becomes 

C(f) = i^)- / 2 !^!/)- 1 fd D x^g~ f ids{isY- 1 - D,2 e- tsi - M2 ~ lt) F{x,x;is). (17) 



J Jo 

The integral over s in the above expression is divergent at v = because of the singular behavior of the integrand as 
s approaches zero. This divergence actually exists for all Re 2/ < D/2. [We will assume throughout this paper that D 
is even] . We therefore regulate the zeta function for these values of v by defining the zeta function for Re v < D/2 as 
the analytic continuation of the zeta function for Re^ > D/2. To this end, we perform D/2 + 1 integrations by parts 
before analytically continuing, to get 

cm = -i{-^r D/2 - ( — — p. — / d D x^- g nidsw x 

{v-D/2)(v-D/2 + l)---v J J 

aD/2+1 . . 

4sj^( e ~ iSiM ' ie)T{X ^ is) )' (18) 

which is regular in a neighborhood of v = 00. This definition of the £ function therefore leads to a regularised one-loop 
effective action. From Eqs. (§) and @, we get 

W 1 = {2(to) D ' 2 (D/2)!(- 1 J d D T^ 9 1 (7 + ln(,, 2 ) + ^ + ^y^y + • • ■ + l) * 

^(-" <M "-' T (— >)„„ 
"Jo (e-^-^F^x;^))!, (19) 

where 7 = Euler's constant = -4^[T(u + 1)] v =q. 

Fixing the renormalization scale /i is equivalent to fixing the (constant) phase of the out-vacuum relative to the 
in-vacuum in flat space, which can be chosen arbitrarily. We will find it convenient to define a rescaled version of fx 
by ln(/i 2 ) = ln(/i 2 ) + -jyj^ + D ^_ 1 + • • • + 1. Then the one-loop effective action becomes 

QD/2 



W 1 = {2{Ak) d ' 2 {D/2)\}^ J dPx^gUy + l n (^))^|-— - ^ m2 ~^F(x, x; is)) 



s=0 



n 



idsHis) 9{is)D/2+1 (e^ ~ ie ^F(x, x; is) J j . (20) 



We may now substitute the -R-summed Schwinger-DeWitt expansion (|15|) in the above formula for the effective 
action and perform the required differentiation and integration term by term. In doing so, we use the integral identity 







00 If 1 1 

idsln(is)(isTe-^ M2 -^ = - (M2 ^. e)p+1 |ln(M 2 - «) + 7 _ 1 - -j , (21) 

for integer values of p. This procedure finally yields 

W 1 = {2{At,) d ' 2 }- 1 j d D x^-g{h{x) + I 2 (x) + I 3 (x)} , (22) 



2 We may equally well analytically continue in D instead of v. 
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where 



oo D/2+1 (_-i\p l+p-D/2-1 

h = -(D/2 + l) E hHM^-ie)-^ £ P U D 2 + 1 - P Y £ n ~" (23) 

l=D/2+l p=0 V ' n=l 

D/2+1 D/2 / l+p-D/2-1 

/ 2 = (D/2 + i)£ £ ( 2 - + ' 7 + MM'-*)- E 

p=l i=D/2+l-p v ' l ' \ n=l 

J t l\{M 2 -ie)- l+D / 2 . (24) 

Z 8 = ( 7 + ln(A 2 )) E L -^/i>/2-p- (25) 
p=0 p ' 

In two spacetime dimensions, the above formulae lead to the one-loop effective action 



1=2 

Similarly, in four spacetime dimensions, we get 



W{ D=2) « {Sir)- 1 J d 2 x^—g^M 2 (in | ^- | —in8(—M 2 ) 

l\ 

(M 2 - ic) 



°° 7! _ 1 

Ew-r'-n ^ ■ (26) 



* = - 2 E A (JP : fe)M + 2 ( Z - - V - 2 )"} ( 2? ) 
h = -p 2 - ~(7 + In | M 2 | -me{~M 2 ))(2] 2 + M 4 ) (28) 
7 3 -i(7 + ln(M 2 ))(M 4 + 27 2 ), (29) 



and the one-loop effective action takes the form 



Wl D=i) « -(647T 2 )- 1 y dS^|(M 4 + 2/ 2 ) (in I ^ I -ztt0(-A/ 2 )) + 3/ 2 

+ E(r 1 + 2(7 - l)- 1 - (l - 2)- 1 ) (M2 _- i£) ,_ 2 /, |> . (30) 



Z=3 



This is the i?-summed form of the effective action. The Gaussian approximation amounts to keeping only the term 
multiplying M 4 in the above series. The terms involving f 2 are required to obtain the correct trace anomaly. 

The first two terms in curly brackets in Eq. ( p0| ) agree with the earlier dimensional regularization results of Parker 
and Toms, except for the step function term 9(—M 2 ) which is purely imaginary and implies particle production when 
M 2 becomes negative (Parker and Toms assumed M 2 > 0). This step function term originates from the identity 
ln(a; — ie) = In | x \ —in8(—x). Within the ii-summed form, therefore, the particle production takes on a very simple 
form, with M 2 = being the threshold for vacuum instability, or the creation of particle pairs. It is conceivable that 
there are physical situations where the imaginary part of the first two terms in the above formula for the effective 
action very closely approximates the actual gravitational particle creation. In this paper, we do not deal with this 
issue further since it does not affect our results. 

Note that Eq. (|3^) is only an asymptotic series expansion in inverse powers of M 2 , arising out of an expansion of 
the heat kernel which ignores, for example, terms that have essential singularities at s — 0. Only the first two terms 
(i.e. up to f 2 ) are necessary for renormalization and for the correct trace anomaly. Also, these terms include the 
convergent infinite sum involving the scalar curvature. Therefore, the approximate effective action based on these 
first two terms is sufficient to indicate the non-perturbative effects coming from the infinite sum of scalar curvature 
terms. This is the form of the effective action that we employ. 
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III. RENORMALIZATION AND OBSERVABLE GRAVITATIONAL COUPLINGS 



Although we are dealing with a free field theory here, the logarithmic dependence of the effective action on the 
curvature in Eq. (^0|) leads to non-trivial effects in strong curvature regions |J . Before we go on to a discussion of such 
effects, it is necessary to understand the contribution of the one-loop effective action to the full gravitational action 
in regions of low curvature. Also, we will now consider the possibility of having multiple particle species contributing 
to the one-loop effective action. For particles of spin | and 1, it has been shown in Ref. [|| that all terms involving 
R in the heat kernel can be summed in a similar manner to the spin case, in a simple exponential form. We will 
therefore consider a generalization of Eq. (pQ) to the form 



W 1 = -U^ 2 )- 1 ^^ / iV=5 {{Mf + Pi f 2i ) In \Mf/p? i \+q i f 2i + ...}, (31) 



where we have inserted arbitrary, species-dependent coefficients rii, pi and qi, and a factor of h is now made explicit. 
The imaginary term is not shown because we are interested here in the part of the effective action that corresponds to 
vacuum polarization. In the low curvature limit considered in this section, we can assume that M 2 — mf + > 0, so 
there is no imaginary term. In later sections where M 2 can be negative but the curvature is very small with respect 
to the Planck scale, the created particles make a negligible contribution to the classical matter already present in 
the Robertson- Walker universe under consideration. The terms of higher order in the curvature correspond to the 
asymptotic series sum in Eq. (|3(]) and do not appear in Eq. (|3lj). For higher spin fields, the heat kernel (and the 
Green's function) generally has a matrix structure. In evaluating the effective action in such cases, one performs an 
additional trace over all internal indices. It will be understood that such a trace has been carried out before arriving 



at Eq. (30). f 2 will be understood to mean the trace of the modified second Schwinger-DeWitt coefficient. The values 
of species-dependent coefficients for spin 1/2 are £ = 1/12, m — —4, pi — 1/2, and qi = 3/2. 

At low scalar curvature <C to 2 ), one may expand the logarithm in Eq. ( |3l|) in powers of R. Noting that f 2i 
can be expressed as the linear combination 

J 2i = a % DR + biR aP R^ + c t R aPj5 R a ^ s , (32) 

the leading terms in the 1-loop effective action then give 

W 1 --ft^TT 2 )- 1 J d i x^g^i\n(my^) + mri l R{l + 2\n(m 2 lfi 2 )) 

+ 1] R 2 (3/2 + ln(m 2 //i 2 )) + {hR^R^ + Cl R aPjS R a ^ s ) ( Pi ln(m 2 //i 2 ) +%) + ...}, (33) 

In the above expression, we have allowed for different renormalization points characterized by the different mass scales 
fii. Changing these mass scales give rise to terms that can be absorbed into the bare gravitational action, as will be 
discussed below. However, one can still remove the /ii-dependence of the full effective action by using our knowledge 
of the observed gravitational coupling constants. To be precise, consider the bare gravitational action 



W g = I d A x^f^g {(k + Sk)(R — 2(A + 6 A)) + {a 1 +8a 1 )R 2 + (a 2 + 5a 2 )R fl „R^ 

+ (a 3 + Sa 3 )R^ 5 R^ s } (34) 



where the counterterms are 



5A = -(US^n^hJ^^t 

i 

5k = (32^ 2 )- 1 ?i^n l TO 2 ^ln( M 2 / / i 2 ) 

i 

fo 1 = (64^ 2 )- 1 ?i^n i e-ln( Ai 2 /M 2 ) 

i 

Sa 2 = (64^ 2 )- 1 a^7i J 6 iP an( A1 2 /^ 2 ) 

i 

6a 3 - (64^ 2 )- 1 ?iJ]n l c l pan( M 2 /^ 2 ). (35) 
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Mi, 2, 3,4, 5 are arbitrary constants of dimension mass, and the /^-dependence of the counterterms is required, following 
||, to cancel the /^-dependence in the 1-loop effective action^]. Adding Eqs. (|33| ) and ( |34|) we obtain the full effective 
action at low curvatures, as 

W = W g + Wx = J d 4 x^=g {-2k A + n a R + a lo R 2 + a 2o R afi R afi + a 3o R afij sR a(3j5 + ■■■}, (36) 

where the subscipt o refers to the observed gravitational constants at low curvatures. These are combinations of the 
bare and induced constants, independent of fli, and given by 

- 2k A = A - ft(647r 2 )" 1 n i m t ln(m?//U?) 

i 

n = k- Ti(32iT 2 y 1 ^2n l m^ i (ln(m 2 /f4) + ~ 

i 

a Xo = a x - n(64 7 r 2 )~ 1 ^ mti (ln(m 2 /n 2 3 ) + ^ 
a 2o = oi2 - ?i(647r 2 ) _1 ^2 n ibi (Pi m ( TO i + q t ) 

i 

«3o = "3 - ^(647r 2 ) _1 ^ n^Cj ln(m 2 //Ug) + <?;) (37) 

Note that one can always absorb the dependence of the observed constants on the fa's into the bare constants. We thus 
have some freedom in shifting terms within the above equations. However, as already stated, we may, in principle, use 
our knowledge of A , k etc. in Eq. ([36]) at low curvature to explore the theory in regions of high curvature. Indeed, 
the full effective action in regions of high curvature now depends only on the observed values of the gravitational 
coupling constants at low curvature, on the physical particle masses rrii and on the values of (these are fixed for 
higher spin fields). No other parameters enter into the effective action. To see this, we add Eqs. (|3l| ) and (|34|), use 
Eq. (|3^) to substitute for the counterterms, and finally use Eq. (|37|) to eliminate the five mass scales in favor of the 
observed constants. In this procedure, the jli dependence in Eq. fl31|) and the \Xi dependence in Eq. ( |37|) cancel the 
corresponding dependences in the counterterms of Eq. ([35]) . Then the full effective action below threshold is given by 

r f y M 2 

W = W g + W x = J d 4 xy^ l-2 Ko A - h—Y,n t mt\n\ — |- | 

Is i 1 

- r bs? E n ^ ln 1 % 1 /2i + ( ai ° + *dU S »&) R2 

+ a 2o R aP R a P + a 3o R af}lS R a ^ 5 + ...}, (38) 



where f 2 i is given by 



Here, 



Pif2i=Pif 2 i + l£R 2 - (39) 



3 In zeta function regularization, the divergent pieces of the one-loop effective action have been thrown away beforehand, 
which is why it is not necessary to introduce those divergences into the counterterms either (although this procedure of 
introducing counterterms could be bypassed in zeta-function regularization, it is necessary in other regularization schemes, 
such as dimensional regularization). In dimensional regularization ||, one explicitly keeps track of the divergent pieces and 
introduces corresponding divergences in the counterterms. Since the divergences are ultimately canceled, one finally ends up 
with Eq. (an) in any case. 
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n = (16ttG )-\ 



(40) 



where G is Newton's gravitational constant, and A D is the usual cosmological constant. 

Eq. (pq ) is a new result, relating behavior at high curvatures to values of the gravitational coupling constants at 
low curvatures. For scalar fields, £ is the only undetermined parameter in the effective action. For higher spins, the 
corresponding parameters have fixed values. We emphasize again that terms involving / 3 and higher, omitted in Eq. 
(|38|), constitute an asymptotic expansion in inverse powers of M 2 and are therefore not necessarily expected to be 
physically significant. The terms retained in Eq. (|3^) are the minimal set of terms necessary for renormalization and 
incorporate the sum of scalar curvature terms in the propagator, ft is readily checked that Eq. ( |36| ) constitutes the 
first few terms in the low-scalar-curvature limit (| £R |<C rn?) of Eq. (|38|). 

In subsequent Sections, we consider the gravitational field equations for a single scalar field, as derived from Eq. 
(|38|), and their consequences. Comparison of Eqs. (|3l]) and ( p0| ) gives, for scalar fields, i%i = l,Pi = 2,qt = 3. In 
Section V, we study the case of constant curvature spacetimes. In later sections, we generalize the investigation to 
Robertson- Walker universes containing matter and radiation. In the next Section, we clarify the usage of the terms 
"perturbative" and "exact" in this paper. 



IV. MEANING OF EXACT AND PERTURBATIVE SOLUTIONS 

In the previous section we derived the effective action, Eq. (B8h, appropriate to the i?-summed propagator. All 



terms multiplied by fi in Eq. (38) will be called "quantum- vacuum" terms, also referred to as one- loop terms, because 
they arise after vacuum fluctuations of the field have been integrated out in a path-integral formulation of the theory. 
However, note that factors of h must also be included in the argument of the logarithmic terms for dimensional 
reasons. Indeed, one has 



In | M 2 /m 2 |= In | (m 2 + fi 2 £R)/m 2 | . (41) 

An approach that is completely perturbative in h would then involve expanding out the logarithmic terms. However, 
the factor of fi 2 appearing in the above equation has its origins in the generalized Klein-Gordon equation 

17l 2 \ 

-a + -^+^RU = 0, (42) 

whose solutions contribute to the tree-level effective action, as opposed to the quantum-vacuum or one-loop effective 
action. That is, the h 2 factor above does not arise from integrating out vacuum fluctuations. Furthermore, the 
quantity m 2 /Ti 2 is the square of the inverse Compton wavelength of the field, and need not be large relative to the 
curvature scalar. Therefore, it is reasonable to regard the logarithmic terms as non-perturbative when expanding in 
h. It is convenient to avoid explicitly inserting factors of fi in the arguments of logarithmic terms (more generally, in 
any term which contains M 2 ), with the understanding that the mass of the field is interpreted as an inverse Compton 
wavelength. 

A perturbative analysis in h 7 as used in Section VI of this paper, therefore treats only the quantum-vacuum cor- 
rections in a perturbative fashion, while keeping the logarithmic terms intact. It is not inconsistent to do so, since 
it defines a regime in which the semiclassical corrections to the effective action are much smaller than the tree-level 
terms involving A„ and k , yet in which the scalar curvature is allowed to be of the same order as the square of 
the inverse Compton wavelength. On the other hand, an exact analysis, as used in Section V below, treats even the 
quantum corrections in an exact fashion. 

The reason we carry out both perturbative and exact analyses is as follows. Exact solutions play a significant role 
when the perturbative analysis breaks down (as signaled by a rapid growth in the contributions of quantum corrections 
to the metric). In order to obtain an understanding of the full dynamics of the metric, it is therefore necessary to 
construct both perturbative and exact solutions to the semiclassical Einstein equations. 

It is well-known [[l5| that semiclassical equations of motion can be perturbatively reduced in a manner such that 
the resulting equations only admit solutions that are perturbative in h. It has been further argued that such solutions 
are the only physically viable solutions of the full semiclassical equations since they do not exhibit runaway behavior 
in the classical limit. 

However, in Section V, we argue for the physical significance of our exact solutions. This argument is based on 
two observations. First, owing to the presence of a mass scale in the theory, h explicitly enters into the semiclassical 
equations only via the dimensionless ratio r = m 2 /m p ,. This ratio is not necessarily small, although perturbation 
theory assumes that it is. Thus, there could be physical effects at large r, not encountered in perturbation theory. 
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Secondly, we will find exact solutions that cannot be expanded in ft and have a well-defined limit as r — * 0, i.e. 
they do not possess runaway behavior in the classical limit. Such solutions must therefore be regarded as physical 
solutions. 

For these reasons, we believe that exact solutions that do not arise from perturbative reduction must be included 
in a complete analysis of solutions of the semiclassical Einstein equations, at least when a mass scale is present in the 
theory. 



V. EXACT VACUUM DE SITTER SOLUTIONS 



In this section, we will consider exact solutions to the equations of motion specialized to de Sitter spacetime. These 
equations simplify considerably for spacetimes of constant curvature. This simplification allows us to include in the 
effective action terms involving / 2 and R 2 . 

As we shall see, there is a rich variety of constant curvature de Sitter solutions, with the scalar curvature being 
highly sensitive to the value of £. We will argue that these solutions are physically viable although they do not appear 
in the perturbative reduction approach. 

Consider, therefore, the effective action of Eq. ([38]), specialized to a single scalar field with mass m and curvature 
coupling £: 

r f i m 2 

W = J d?xV=9 |-2k»A„ - h— 2 m* In | -j | 

• + h it^ (> - 21 ° i f£ i)) r - R ^ ta i £ i a + <w< 2 « 2 } ■ («> 

where we have assumed that the constants 0t\ o , 0L2 O and a^o are negligibly smalQ 

One may now use Eqs. (B1-B9) of Appendix B to obtain the equations of motion resulting from the variation of 
Eq. (fl3|). When specialized to a constant curvature maximally symmetric spacetime, in which 

R 

R^ap = -^{g^agv/i - 9^ a g^fi), (44) 
these equations yield a single algebraic equation satisfied by the scalar curvature 

2M4Ao - R) _ - JL {,„. ta 1 £ 1 (, + ? JL) _ l„ m (, + «£ 

-f^-h^-Jf-^)]. (45, 



M 2 2^M 2 V 1080, 

Note that we recover Starobinsky inflation by first taking A Q — > and m — ► 0, and then taking the lim it f —> 0. 
The resulting solution has scalar curvature R = 69120 tt 2 k TtT 1 , in agreement with Starobinsky's results |lq| . The 
Starobinsky solution is non-analytic in ft, and is not well-defined in the limit Ti — > 0. The existence of such solutions 
has motivated arguments in favor of the perturbative reduction scheme ]lq] , which discards such solutions in a self- 
consistent manner. 



4 For the purposes of this section, it is actually sufficient to assume that terms involving at a , 0120 and oizo combine to yield 
ao/2+negligible contributions, where a is some constant, because the variation of J d i x y f—gf2 vanishes in a constant curvature 
spacetime. 

5 The right-hand-side (RHS) of Eq. ( J45| ) gives the trace of a stress-tensor in de Sitter space which is not the same as the 
stress-tensor for the Bunch-Davies vacuum state. The stress-tensor on the RHS of Eq. (^) corresponds to a different state and 
is defined by variation of W . The leading terms of the stress-tensor agree with those arising from the Gaussian approximation 
of Bekenstein and Parker which is known to be a good approximation in de Sitter space, particularly for closely spaced 
points as are used in defining the stress tensor. The additional term, (1/2)(1/1080)£(-R 3 /M 2 ) in Eq. (Uft) yields the correct 
trace anomaly. Furthermore, this stress-tensor is conserved, and vanishes in flat spacetime (R = 0). The RHS of Eq. Jisl ) 
diverges at M 2 — m 2 + £R = 0. Similar divergences also occur in the stress tensor in the Bunch-Davies vacuum state in de 
Sitter space, which diverges at m 2 + (£ + n(n + 3)/l2)R — 0, where n is a non-negative integer. 
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However, such arguments weaken in the presence of an additional mass scale m in the theory. To see this, consider 
Eq. ([l5| ) rewritten in terms of the dimensionless variables 

y=4 ( 46 ) 

ft- 1712 r,~s 

r = 77i ■ ( 4? ) 

lOTTKo 

Recall that m refers to the inverse Compton wavelength of the particle, so y and r are indeed dimensionless. In terms 
of the actual mass of the particle, r is given by 



actual 



'PI 



(48) 



where mpi = \/1Qttk Tl. 

Eq. (ft5| ) then takes the form 



(1 + &) _ 1 1 + & 1 {1 + \iy + ly (r - jij) } - £ (, - £ j . («) 

The solution for ?/ is a function of three dimensionless parameters: r, £ and A /m 2 . In a perturbative reduction 
approach, r is regarded as a small parameter and the solution for y is constrained to be an analytic function of r. 
However, it is plausible that the limit Ti — > does not imply r — > 0, i.e. that the mass m is rescaled in such a manner 
that r stays constant (or roughly constant) as % — > 0. For example, in string theory, one would expect that all masses 
(including the Planck mass mpi) are generated by a single scale (the string scale), in which case the dimensionless 
quantity r would be independent of h. In the absence of knowledge of the precise nature of a fundamental unified 
theory, it is therefore prudent to consider the possibility that the parameter r is some finite quantity, not necessarily 
small, and to treat it in a non-perturbative fashion^. In subsection B, we give a further justification for the physical 
validity of solutions that arise without expanding in r. Namely, we find that for £ < 0, these solutions have a 
well-defined classical limit as % (or r) — > 0. 

We now consider numerical solutions to the algebraic equation ([l9]). This equation, in general, has solutions with 
both positive and negative scalar curvature. We will, however, focus on the positive scalar curvature solutions (y > 0), 
corresponding to an inflating de Sitter universe. Figs. 1 and 2 are plots of the left hand side (LHS) and right hand side 



(RHS) of Equation (49), as functions of y, for various values of the three parameters r, £ and A Q /m 2 . The point(s) 
of intersection of the LHS and RHS correspond to solutions for y. These plots are convenient ways of identifying the 
solution space because the LHS depends solely on the parameter £, while the RHS is a linear function of y , with slope 
given by r, and intercept given by A /m 2 . In all plots, therefore, the straight line is the RHS of Eq. (f4£fr). Increasing 
the value of r will decrease the slope of the straight line. Increasing A D will shift the straight line up. It is convenient 
to consider two ranges of values of £ that give qualitatively different behavior of the LHS of Eq. (|4^). These are: a) 
£ >0, and b) I < 0. 

We will find that there exist solutions with non-zero scalar curvature, even if A D = 0. For £ > 0, the most interesting 
solutions of this type exist for values of £ very close to (1080)~2 . For other values of £ > 0, there are either no solutions 
or solutions for which the scalar curvature is of order m Pl , which may thus be unphysical. 

For £ < 0, we find that there exist solutions with R ~ — m 2 /£, for a large range of values of £ and r, and for small 
and large values of A D /m 2 . These solutions are of greater interest for the purposes of this paper, and the reader may 
safely skip directly to subsection B. 

A. Solutions with £ > 

For £ > 1080 -1 / 2 and A Q = 0, the only solution to Eq. ( f49| ) is the trivial solution y = 0, because, for y > 0, the 
LHS is always negative, while the RHS is always positive. However, addition of a non-zero value of A Q allows the RHS 



6 DifTerent arguments for the physical significance of solutions that do not arise from perturbative reduction have been given 
by Wai-Mo Suen O. 
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to take negative values, and leads to a non-trivial solution with a value of y slightly lower than the classical value 
4A Q /m 2 . The deviation from classical behavior will typically be very small due to the extreme flatness of the LHS 
graph near the origin. 

As one lowers the value of £ towards (1080) -1 / 2 , the LHS acquires a local maximum (Fig. [|). The position of this 
maximum is very sensitive to the value of £ when £ is close to the value (lOSO)^ 1 / 2 . For large values of r (r = 10 in 
the figure) , a non-trivial intersection with the RHS graph is now possible, with y- values ranging from about a hundred 
to extremely large values, depending on the precise values of £ and r. The presence of a non-zero A Q term would shift 
the RHS graph down and increase the value of y even further. 

When I = (1080)- 1 / 2 = 0.03042 . . ., the non-trivial intersection point of the two graphs will occur at an extremely 
large value of y, for typical values of r. One may obtain an analytical estimate by considering an approximation to 
Eq. ( fi"9| ) for large y, after setting £ = (1080)~3. This gives 

]xi-^= = —Vvm+l. (50) 
VT080 r 2 

For r — 1 (to = mpi), the above equation implies a value of y approximately equal to about 10 91 . This value increases 
as r decreases (i.e. by decreasing m and holding mpi constant) and vice- versa. However R = m 2 y = m^ry will acquire 
a minimum value for some r > 1. Also, addition of a A Q term tends to decrease the value of y and, consequently, R. 

For <J < (1080)- 1 / 2 , there is a non-trivial solution with an extremely large value of y. An analytic approximation 
to Eq. ( [49|) can be made once more, after noting that the term involving y 2 in the LHS is now the dominant term 
(this term vanishes when £ = (1080)~ 1//2 ). Thus we obtain, in this regime, with zero cosmological constant, 

'-TljSso-O ■ (51) 

y therefore scales linearly with r _1 . However, the scalar curvature itself is essentially independent of to in this regime, 
and is given by 

R = 4 ™"(Tk6- f ) 1 (52) 

The two equations above predict that y (or R) will decrease as £ decreases. This behavior is borne out by numerical 
study. However, when £ is extremely small, £?/ can become small, and the large-£y approximation breaks down. For 
very small values of £, the solution for y then increases with decreasing £, giving y — > 00 as £ — > + , as expected from 
Eq. ©. 

To summarize, for £ > and A Q = 0, non-zero solutions for constant scalar curvature occur for £ values very close to 
(1080)" 1 / 2 (this corresponds to the curvature coupling £ being very close to the conformal fixed point £ = 1/6, because 
(1080)^2 is a small number). Such solutions also occur for £ = (1080) _1 ^ 2 ; however, in this case, the curvature is 
typically many orders of magnitude greater than Planck size, and may thus be unphysical. 



B. Solutions with £ < 

For £ < 0, the LHS of Eq. fl4S| ) becomes singular at y = — £ . This value of y becomes an exact solution to Eq. 

( ft9| ) in the "classical" limit r — > 0, as can be seen by letting y = — £ + e(r) in Eq. (Q), and showing that e(r) — > 
as r — > 0. This fact constitutes another argument against the perturbative reduction scheme in this case, because the 
scalar curvature does not exhibit runaway behavior in the classical limit. As we shall see, most solutions with r < 1, 
will lie very close to the value y = —t; . Significant deviations from this value will occur only for very large values 
of r, or when a non-zero A D term is present. 

Fig. H is a plot of the LHS and RHS of Eq. ( |49| ) with representative values £ = —0.03 and A Q = 0. An exaggerated 
value of to = 4top; is used to fully display the graph. However, as we will see, the solution for y is largely insensitive 
to the precise values of £ and r. 

The straight line LHS graph intersects the RHS graph at a y-value very slightly larger than — £ (= 33.3 ... in this 
case). The corresponding scalar curvature is given by R ~ to 2 £ . The steep vertical ascent of the LHS graph near 
the value y = — £ implies that the solution for y is not very sensitive to the precise value of r (note that decreasing 
r increases the slope of the straight line, giving a solution even closer to the value £ 1 ). For most values of r, except 
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for exceptionally large values, we have a solution at R = — m 2 £ plus a small correction. The presence of a A Q term 
has the effect of shifting down the RHS graph, leading to larger values of the scalar curvature. However, in this case, 

a second solution will appear with y < — £ , the y- value being very close to the classically expected value 4A Q /m 2 , 
for a small cosmological constant. As the cosmological constant is increased to large values this second solution will 
now approach y = — £ , while the first solution now corresponds to the classically expected large value of y. The 

solution y ~ — £ is therefore fairly robust; it exists for large ranges of the parameter r and for both small and large 
values of A D . Numerical investigation shows that such a solution also exists for a large range of values of £. 

It should be noted that a y value slightly larger than — £ corresponds to a small negative value of M 2 , and 
therefore the effective action evaluated on such a solution acquires a small imaginary part giving rise to a small rate 
of particle production. In the regime we consider in this paper, the effects of this particle production are negligible. 

The qualitative behavior of the LHS of Eq. ( |49| ) for £ = — (1080) - ^ is very similar to the behavior for £ = — 1, 
and quite insensitive to the precise value of r, as in the previous case. The solution will now be very close to 

y = — £ ~ 32.86 .. ., implying that the scalar curvature is about an order of magnitude larger than m 2 , for a large 
range of values of m. Again, this is a robust solution as discussed earlier. 

To summarize, for £ < 0, we find physically reasonable values (in the sense that they could be small with respect 
to m 2 Pl ) of the scalar curvature, approximately given by 

o 

R~ =r- = to 2 (53) 

for a large range of values of r, £ and A D . The approximation in Eq. (^3|) breaks down if A Q = and r 3> 1, in which 
case the only solution corresponds to an extremely large scalar curvature (i? 3> m 2 Pl ) whose precise value depends 
on r and £. This large R solution will give rise to a large imaginary contribution to the effective action, and so one 
expects copious amounts of particle production to occur, which could, in turn, bring the scalar curvature down to 
reasonable values. This would be consistent with a gravitational Lenz's law mechanism 
this issue here. Instead, we now turn to a perturbative analysis of the semiclassical Einstein equations. 

VI. THE GROWTH OF QUANTUM CORRECTIONS TO THE SCALAR CURVATURE 

In this section, we will analyze, the effect of logarithmic curvature terms in the effective action on a Robertson- 
Walker cosmology. We analyze this effect in two ways, valid for spatially open, closed and flat models. First, we 
consider a universe with mixed matter and radiation and assume that the scalar curvature is slowly varying so that 
its derivatives can be ignored. This analysis shows that perturbative (as defined in Sec. IV) quantum corrections to 
the scalar curvature can sharply increase in magnitude near a time tj that is determined by m and £. In particular, 
the existence of an ultralight mass in the theory can lead to significant quantum effects close to the present time. 
Since the background classical scalar curvature is decreasing, the effect of quantum corrections is to prevent the scalar 
curvature from decreasing further. However, as quantum corrections become more and more significant, perturbation 
theory breaks down and cannot be relied upon to give the full behavior of the scalar curvature. We then carry out a 
second analysis of the behavior of the scalar curvature for all times after tj , without using perturbation theory. This 
second analysis indeed reveals that, for t > tj, the scalar curvature tends to decrease extremely slowly, ultimately 
approaching a constant value. The analysis thus displays consistency with the original assumption of slowly varying 
scalar curvature. 

In Section VII, we will use the behavior of the scalar curvature to construct a model universe in which a matter- 
dominated cosmology transits to a mildly inflating de Sitter cosmology at the time tj. 

A. Quantum Corrections to the scalar curvature 

In this subsection, we will consider a classical Robertson- Walker cosmology with mixed matter and radiation, and 
treat the quantum effects involving logarithmic curvature terms in a perturbative fashion (recall the discussion in 
Section IV, which shows why an expansion of the logarithm itself is not appropriate in a perturbative treatment 
of quantum- vacuum terms). The universe will deviate from classical behavior when the quantum effects become 
sufficiently large. The essential idea is to allow for the possibility of quantum effects being significant at the present 
time. We will find that this is possible if there exist very light mass fields with £ < 0. 



18 HJ]. We do not address 
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Our starting point for the analysis here is the effective action of Eq. (pig)). Variation of this effective action, 
specialized to a Robertson- Walker metric, will generally yield terms involving derivatives of the scalar curvature R, 
in addition to the terms of Eq. p5|). In this subsection, we assume that terms involving derivatives of the scalar 
curvature are negligible. We will find that, for light mass fields (m -C mpi), this assumption is justified because 
derivatives of the scalar curvature remain small until the magnitude of the quantum-vacuum contribution to the 
scalar curvature itself becomes comparable to the classical contribution to the scalar curvature. Beyond this point, 
the perturbative analysis breaks down. In the next subsection, we carry out an exact analysis, valid for t > tj. 

Although we ignore derivatives of the scalar curvature when carrying out the variation of Eq. (|43|), the resulting semi- 
classical Einstein equations differ from Eq. ([l5]) in two important respects. First, the fact that the Robertson- Walker 
universe has fewer symmetries than de Sitter space gives additional terms in the semiclassical Einstein equations. 
Secondly, we will include a classical stress tensor source representing mixed matter and radiation. We also set A G = 0. 
The trace of the semiclassical Einstein equations, expressed in terms of dimensionless variables, then takes the form 
of a simple generalization of Eq. ( pif ) : 

1\ , 3 - , 1 .ss (t* 1 



r(i + iy) in 1 1 + e?y i - r— *a=- i + + k 



l + £ y { 2"" 2" V 1080 
T 
2m 2 K 



^[v+T^r), (54) 



where T is the trace of the classical stress tensor, and v is a quantity which vanishes in de Sitter space: 

' ( -R 2 - R^R^) . (55) 



~ 180m 4 \A ^ 

All quantum contributions are grouped in the LHS of Eq. (|54|), while its RHS contains classical terms. 

Consider now the full semiclassical Einstein equations with a classical stress tensor source representing mixed matter 
and radiation, 

G>iV = i ( ( Pm + t Pr ) """" + l Pr9 ^) + ° {h) ' (56) 

where p m and p r represent the matter and radiation energy densities, respectively, and 0(K) represents the quantum 
contributions. The equation above implies 



R_ _ Pm Rq 

to 2 2n m 2 m? 
(p m + (4/3)p r ) 2 



— = — 2+2' ( 57 ) 



960m 4 K 2 



0(h), (58) 



where Rq, as defined by Eq. (p7|), is of order h. In a treatment perturbative in h, we replace y and v in the LHS 
of Eq. (pj) by their classical values using the above equations, because the factor r is of order K. This treatment is 
valid as long as the LHS of Eq. ([34j) is small with respect to the term irT / (m 2 k ) in the RHS. In the RHS of Eq. j5^ ) 
we keep the quantum- vacuum contribution to y as well (i.e. the term Rq/iti 2 ). Furthermore, in a Robertson- Walker 
cosmology with metric 

ds 2 = -dt 2 + a(t) 2 ( + r 2 dn 2 ) , (59) 
\ 1 — kr z ) 

where k = —1, +1 and denote spatially open, closed and flat universes respectively, we have p m oc a~ 3 and p r cx a~ 4 . 
It is then more convenient to express y and v in terms of the present matter and radiation densities p mo and p ro , and 
the redshift z. We therefore introduce new dimensionless variables d m and d r , given by 

e - dm = (60) 
z- dr = Z^-' (61) 



and the redshift z given by 
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(62) 



where clq is the scale factor at the present time. Furthermore, owing to the fact that we will study the effect of very 
light masses, we also redefine the ratio m 2 /m|> ; : 

r = e- s . (63) 



Eq. (|54|) can now be rewritten by substituting for y and v using Eqs. (]57]) and (|5S|) with the redefined variables of 
Eqs. (|60D , (|6C|), ( |62| ) and ( |63| ) above. The resulting equation, correct to leading order in h, gives an expression for the 
quantum correction to the scalar curvature, Rq, as a function of redshift. Dividing by R c i — p m /(2fc ), we obtain 

Rq _ Rq 



Rcl Pm/(^K ) 

1 ^(l + z) 3 e s - d ™ ( , 3- n , , 3 , 1 ^ 2 1 



I + + zfe s - d -j In | 1 + -£(1 + z) 3 e s - rf >" | 

|l + |f (1 + Z )V^ + 1 (j 2 _ _L) (1 + ,)6 e 2(5-, m ) 



2 1 + + z)3 e s~<i. 

_ J_(l + ^7 e 25-d r -d m _ J_/j + z) 8 e 2(S-d,.)U 

480 v ; 960 v ; J J 



(64) 



We wish to find the redshift ranges for which the quantum contribution to the scalar curvature is significant in 
comparison with the classical contribution (i.e. Rq/R c i of order 1). The values of d r and d m in Eq. d&jj ) are 
determined by the radiation and matter energy densities at the present time. Black body radiation at a temperature 
of 2.726 K gives a radiation energy density p r0 ~ 7.81 x 10~ 34 g/cm 3 f2(||I|. This gives a value for d r ~ 288.06. The 
matter density is not known to good precision. With the conservative estimates Q > 0.1 for the ratio of the present 
matter density to the critical density, and Ho > 50km/ (sMpc) for the Hubble constant at the present time, we obtain 
PmO > 4-70 x 10~ 31 g/cm 3 for the matter density. This gives d m < 281.66 for the exponent of Eq. (|60|). 

For very light masses (S > 150), numerical investigations of the ratio Rq/R c i as a function of redshift z, using Eq. 
(U) with | £ |~ 1, reveal that there are two distinct regimes for which this ratio is close to or larger than 1. The first 
such regime occurs at extremely high redshifts [z > 10 26 ), close to the GUT scale. [In the standard cosmology, the 
Planck era occurs at a redshift of about 10 31 , while the GUT era sets in at a redshift of about 10 26 .] It is expected 
that quantum effects would play a significant role at such high redshifts. The second, unexpected regime for which 
the quantum-vacuum terms become large occurs at relatively low redshifts. This second regime exists only for £ < 0, 



and occurs when the factor 1 + (1/2)£(1 + z) 3 e s dm in Eq. (64) approaches zero. This corresponds to values of z 
near a redshift Zj given by 

(1 + Zj ) 3 = -2r 1 e dm - s . (65) 

In between these early and late regimes, the quantum contribution to the scalar curvature is extremely small and 
slowly varying, and the evolution of the universe is well-approximated by its classical evolution. 

Numerical investigation of the behavior of the ratio Rq/R c i for values £ = — 1 and d m ~ 280, and for very light 
masses (S > 150), further reveals that in the late regime of significance of quantum vacuum terms (i.e. near z = Zj), 
the scalar curvature sharply increases as z — > Zj. Since Zj depends on S (from Eq. (|65|)), and therefore on m, the 
value of to dictates the value of Zj at which quantum vacuum terms can become significant. As we shall see later, 
an ultralight mass can lead to quantum vacuum effects becoming significant at roughly half the age of the universe. 
However, it is important to note that when the scalar curvature begins its rapid increase near redshift Zj, quantum 
effects begin to dominate and the perturbative analysis itself breaks down. 

For our purposes, we will only need the result that the ratio Rq/R c i tends to rapidly increase as z — > zf. The total 
scalar curvature R, which can be written as 

R = Rd + (66) 



R, 



cl 



is thus a product of a decreasing function (R c i) and an increasing function. For z 3> Zj, R ~ R c i decreases, and as 



z 



the quantity Rq / R c i tends to suppress the further decrease of R. 



It will be useful to have an estimate of zj relative to the redshift at matter-radiation equality, z eq . The quantity 
z eq is defined by the condition of equality of matter and radiation energy densities: 
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p m o(l + z eq f = p r0 (l + z eq f. (67) 
Combining Eqs. ©, @, © and ©, we get 



1 + * /2N 1/3 



1 + z eq \ -£ 



l_ — I _ | 4d m /3-d r -S/3 



e «+*/*-*r-*/\ (68) 



As earlier, we suppose that d r ~ 288.06 and c? m < 281.66. Also, we will additionally assume that S 1 > 278.[] With 
these ranges, we obtain 



1 + Zi ( 2 



l + z e 



1/3 



< 5.61 x 1CT 3 . (69) 



Thus, if — £ is of order 1, we have Zj <C z eq . 

It follows, for this range of masses, that the quantum corrections become significant at some time during the matter- 
dominated stage of the expansion. In Section VII, we derive a formula for the cosmic time tj corresponding to the 
redshift Zj, for a spatially open matter-dominated universe. 

To find the scalar curvature at tj, we use Eqs. (|oJ) and (|6l| ) to write 

(1 + 2j )V-. s(1 + 2j )3^ 
Pmj 



= (70) 

where p m j is the matter density at time tj, and the last equality in the above equation follows from the classical 
Einstein equations. Comparing the above equation with Eq. ( |65j ) shows that 

R cl {tj)=m 2 /{-l)=m 2 . (71) 

We will now carry out a second analysis, valid outside the perturbative regime, which shows that the scalar curvature 
indeed approaches a constant value as z — ► Zj (or t — ► tj). 

B. Behavior of the scalar curvature for t > tj 

As t approaches tj , the arguments of the previous subsection show that the classical scalar curvature of the matter- 
dominated universe approaches the value Rj — to 2 . We will now argue that the scalar curvature does not decrease 
further for t > tj, i.e. it saturates to a value very close to, but slightly larger than Rj. That is, we will show that there 
exist approximately de Sitter-like solutions to Eq. (^i|) for which the scalar curvature does not change significantly, 
although the matter density keeps decreasing, eventually approaching zero at very late times. We hnd that, for matter 
densities less than p m j there exist de Sitter type solutions of Eq. (p4) of the form 

R = m 2 {l-e) (72) 

such that | e |<C 1, and R is an extremely slowly varying function of the matter density. We will assume throughout 
that the mass m is very light (r « 1), and that — £ is positive and of order 1. 

To show that Eq. (^) is indeed a solution, we substitute Eq. ( J72| ) in Eq. (|54|) and assume v <C 1, so that one can 
effectively ignore v. Then Eq. (^||) takes the form 

eln | e | - ^ jl+ | + \{l - e) 2 (l - (1080f )~ l ) + 



2vr ( e - 1 T 
r 



£ 2to 2 k 



(73) 



7 This assumption will be justified in the next section, where we find the value of the mass (and therefore S). 
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Define 



T V=-^h-V- (74) 



2m 2 n 2m 2 K, 

-r-1 



If Pm = Pmj, where p m j is found from Eqs. ([35]) and (|7C|), 6 = 0; and if p m = 0, S = — £ . We find that | e | is 
always small within this range of 5 values. In fact, we show that it is small for any positive value of 5. In the small 
| e | approximation, Eq. (|73|) becomes 



-i 



1 \ 2tt 



c 



2(1080)£ 2 



v\= — l= + 5\. (75) 



To estimate v, we assume that it takes a value in between its value for a classical matter-dominated universe at time 
tj, and a de Sitter universe, for which v = 0. In a classical matter-dominated universe, v is given by Eq.(|5^) as 

/4 / 76 n 

At time £j, p m = p m j = 2n rn 2 . Therefore, at tj, we obtain 

Vj = -(240FT 1 . (77) 
We therefore assume that, for t > tj, v lies in the range 



(240CV 1 < « < 0. (78) 



It is convenient to define an additional quantity 



P=±-[v 3 ■ (79) 

2lT V 2(1080)f/ 

For r«l, and v given by Eq. (|78|), it follows that | (3 |<C 1- One may now solve Eq. (]7^) for e. This yields 

e = -ie(5±(^ 2 +4/5r 1 )). (80) 

In order to choose the correct sign in the above equation, we will require that the scalar curvature approach its 
classical value p m /(2K ) for large values of the matter density, i.e. for p m /(2m 2 K ) S> 1. For large matter densities, <5 
is a large negative number. To get the correct value of the scalar curvature, one must then choose the solution with 
the minus sign in the above equation. Even though this value of e is not small, continuity demands that we keep the 
solution with the same sign for all 6. Thus we have 

c = (t - (S 2 - WV)) ■ (81) 



It is clear, for S > 0, that | e | is always a small number, reaching a maximum value of y /3£ at S — (we have assumed 
that r < 1). 

The approximation | e |<C 1, which was made in order to arrive at the solution (|Sl|), is thus valid for t >tj. Within 
this range of time, the value of e evolves from 

efa) = (82) 

when p m = p m] (<5 = 0), to 

e(oo) a -0£ (83) 
when p m — (S — — £ ). The fractional change in the scalar curvature during this range of time is given by Eq. 



17 



because | e |<C 1 during the entire time range. Thus, with Ae = e(oo) — e(tj), we obtain 

| ^ |~ y/ft (l + y/(%\ « 1, (85) 

because r< 1. 

To summarize, we have argued that, even in the presence of matter, there exist de Sitter type solutions in which 
the scalar curvature is very close to the value to 2 , and is very slowly varying, as long as the matter density is such 
that S > 0. For large negative values of 6, i.e. large p m in Eq. (|74|), it is clear that the scalar curvature continuously 
changes to the scalar curvature of a matter-dominated universe. 

Coupled with the findings of the previous subsection, in which we showed that quantum effects near tj tend to 
prevent the scalar curvature from decreasing further, these results strongly point to a cosmology in which a matter- 
dominated universe transits to a de Sitter-type universe. In the next subsection, we outline such a model. 



VII. MATTER DOMINATED EXPANSION LEADING INTO ACCELERATED EXPANSION 

The damping of scalar curvature, which begins at a time close to tj, supports the idea that the universe undergoes 
a transition from a matter-dominated phase to a mildly inflationary de Sitter phase. The arguments of the previous 
subsection show that the de Sitter phase is, to good approximation, described by a solution with £ < and R ~ to 2 , 
of the type found in the Section V, and that the presence of matter does not significantly change such a solution. 
Within a rigorous framework, the matter-dominated and de Sitter phases must be joined in a sufficiently smooth 
manner to gaurantee regularity of all curvature components at the joining point. However, the perturbative analysis 
of the previous subsection shows that the quantum effects become significant over a short timescale, after which the 
scalar curvature approaches a constant value. This effect allows us to consider, for comparison with observation, 
an approximation in which an exact classical matter-dominated solution is joined (at time tj) to a de Sitter solution 
generated by quantum effects. Assuming a spatially open cosmology (k = — 1 in Eq. (|59|)), such a model is represented 
by the following scale factor: 

a {t) — c o sinh 2 (?/>/2), t — —co(smhtp — tp), t<tj 

a(t) = a~ 1 smh(a(t + c 1 )), t>t d . (86) 

Here, tp parametrizes a(t) and t during the matter-dominated stage. For t > tj, including the present time to, the 
universe is in a de Sitter phase. 

Of the five parameters, cq, a, tj, c\ and the present cosmic time to, that characterize the model based on Eq. (pq), 
not all are independent. The scalar curvature at the time of joining, Rj, must be equal to the constant scalar curvature 

during the later, inflationary phase, and is determined by the single scale to = m/ \J — £. This requirement and the 
requirement of continuity of the scale factor at the joining point, constitute two constraints on the five parameters. 
The remaining three parameters can, for convenience, be taken to be (i) the present Hubble constant Hq, (ii) the 
present ratio of matter density to critical density FIq, and (iii) the mass scale to. Here, we express the five parameters 
defining the model of Eq. ( p6|) in terms of these three basic parameters. 

First, the scalar curvature during the later, de Sitter phase is given by 12a 2 . Setting this equal to to 2 , as required 
by the de Sitter solutions of Section V, we get the relation 

(87) 



The scalar curvature during the matter-dominated phase is given by R = Scqci 3 . The classical Einstein equations, 
which hold during the matter-dominated phase, thus imply 

3c - (8TrG)p mj a 3 j = {8TrG)p m0 al, (88) 

where dj and ao represent the scale factor at tj and the present time tg, respectively, and p m j and p m o are the 
corresponding matter densities. The second equality in the above equation follows from the fact that /9 m a 3 is constant 
during the evolution, a consequence of conservation of the matter stress-energy. 
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Thus, Eq. ( pq ) gives 

c /a 3 = M 2 , (89) 

where Qq — 8nG p rrl Q / (3H'q) is the present ratio of matter density to critical density. 
The Hubble constant at the present time is given by Eq. ( p6| ) as 

Hq = a coth(a(to + ci)) 

= a^ l s Jl + ala\ (90) 

where a is the scale factor at the present time t . We solve for a , and use Eq. ( |S7j ) to get 

a = (ff 2 -m 2 /12)- 1 / 2 . (91) 

Combining Eqs. (|9|) and (|9l]), we obtain 

c - n H$(H<- - m 2 /12)~ 3 / 2 . (92) 

We may obtain the scale factor at time tj , aj , by requiring that the scalar curvature of the matter-dominated phase 
approach the value m 2 as t — > t J . This condition yields 

= US- (93) 

Substituting for Co from Eq. we then have 

a, = (3r»oi? 2 M 2 ) 1/3 (^o - ™7l2) -1/2 . (94) 
To obtain tj, we use the matter-dominated solution in Eq. (|S6|) to get 



^,•=2 sinh y o,jC (95) 

and 



tj = ic (sinh^ - Vj) 



c ^y/(a j Co 1 )(l + a j c 1 ) - sinh 1 \[a~c ~^j 



(96) 



To obtain ci, we use the de Sitter solution in Eq. ( |86| ) to get 

ci = a -1 sinh - (era?) — ij. (97) 

Finally, to obtain the present cosmic time ioj w ^ again use the de Sitter solution, which yields 

to = a -1 (sinh -1 (aao) — sinh -1 (araj)) -I- tj. (98) 

All parameters in the model are now expressed in terms of m, Hq and Qq. 

In the next section, we will compare the predictions of this model to recent data from high-redshift Type la 
supernovae B, using magnitude-redshift curves obtained from our model. 



VIII. COMPARISON OF THEORY WITH HIGH-Z SUPERNOVAE DATA; PREDICTION OF PARTICLE 

WITH MASS ~ 1(T 33 EV 

Recent observations of Type la supernovae at high redshifts indicate a negative value of the deceleration parameter 
at the present time, i.e. an accelerating universe ||. Previous attempts to account for this phenomenon invoke a 
cosmological constant, or a classical scalar field, quintessence [^3|, with unusual potentials. To explain the observed 
acceleration effect by means of a cosmological constant, it must contribute a term to the Einstein equations that is 
of the same order of magnitude as that attributed by the present matter density. On the other hand, quintessence 
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models that account for the acceleration effect typically involve potentials that would give rise to nonrenormalizable 
quantum field theories. 

In the model we present here, the value of Hq is fixed by low-redshift measurements, while the remaining free 
parameters m and Qq are determined by the SNe-Ia data. Once the mass scale m is determined, no fundamental 
parameters in the effective action need be chosen to fit the supernovae data. Furthermore, the theory we work with 
arises out of a renormalized effective action. 

Comparison of our model to SNe-Ia data is achieved by fitting calculated magnitude-redshift curves to the data. 
The difference between the apparent magnitude (m) and absolute magnitude (A/) of a source is given in terms of the 
luminosity distance d^ to the source, by 

m-M = 51og 10 ^- + 25. (99) 

The luminosity distance itself is given by 

d L = (l + z)a r 1 , (100) 

where ao is the present scale factor, and r\ is the comoving coordinate distance from a source at redshift z to a 
detector at redshift 0. For Robertson- Walker universes, r% is given by the equation 

dr i f z dz' 

o (i^w = a ° Jo m*v (101) 

where k = 0, +1,— 1 correspond to flat, closed and open universes, respectively. For a spatially open universe, the 
above equation reduces to 

( z dv 1 

sh± "-^L iwr (102) 

Consider a universe represented by the model of Eq. (|S6|), which is a spatially open matter-dominated universe 
prior to time tj , and transits into a de Sitter universe at tj . Let Zj be the redshift at time tj . For z < Zj , the universe 
is in a de Sitter phase, with Hubble constant given by 

H = a coth a{t + c\) 



(103) 



Substituting the above into Eq. (102), and performing the integration, we obtain, for z < Zj, 

ri<w = ^fvTTM-X7(^f), (104, 



where n<(z) denotes ri(z) for z < Zj. For z > Zj, the RHS of Eq. (1103) separates into two contributions: 



«>> = °^f <105) 



where r\>{z) denotes ri(z) for z > Zj. 

During the matter-dominated phase, the Hubble constant is calculated as 



H = \J coa 3 + a 2 
1 + z 



a 



yjl + coa^il + z). (106) 



Using Eq. (103) for z' < Zj and Eq. ( JlOq ) for z' > Zj, the integrations in Eq. (105) may be performed to yield 

^ M -^(^ <M)+1 .( g|;j;;»w^;» )), im 
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where 

g(z) = y/l + oa^il + z). (108) 
Eq. ( |100| ) gives the luminosity distance, d,Li, for this model, as 

dm = (1 + z)a ri < (z), z < Zj 

= (1 + z)a Q r 1> (z), z>Zj. (109) 

For comparison, the luminosity-distance-redshift relation for a matter-dominated Robertson- Walker universe with 
zero cosmological constant and ratio Slo of present matter density to critical density, is 

d L2 (n , z ) = 2# - 1 ft 2 (n oZ + (n - 2)^1 + n QZ - 1)) , (110) 

and for a spatially flat matter-dominated (Qq — 1) universe, 

d L2 (l,z) = 2H Q 1 VlTz~(VTTz~- 1). (Ill) 
It is convenient (see Eq. (^)) to define 

A (m -M)(,) = 5.„ gl0 ( 5 -|a_). (112) 

Fig. H is a plot of A(m — M) vs. z, along with a plot of SNe-Ia data acquired from Ref. ||. The two solid curves 
represent di,{z) = di,\(z). In plotting this quantity, all parameters appearing in dm have been expressed in terms of 
the three basic parameters to, Hq and f2o, using the relations derived in the previous section. Also, the value of Hq 
has been set to 65km/ (s Mpc). Thus there are two quantities, to and Qq, which parametrize the solid curves. The two 
curves shown in the figure give a reasonable fit to the data, and correspond to a) to = 3.7 x 10 -33 eV and flo = 0.4 
(higher solid curve), and b) m = 3.2 x 10 -33 eV and flo = 0.3 (lower solid curve). 

The general features of a family of curves parametrized by (m, Qq) are as follows. For a fixed value of m, decreasing 
f2o has the sole effect of increasing the redshift at which the transition occurs, i.e. a smaller value of f^o will move 
the transition further from the present time. Thus we cannot rule out the possibility that, with more observations at 
higher redshift, a better fit to the data could be obtained with lower values of Hq. However, the data do not allow 
the joining points in Fig. (||) to occur at smaller Zj, so the values of fio shown in the plot do represent a rough upper 
bound on f2 in our model, and lead to the conclusion stated earlier, f2 < 0.4. As is well known, an early inflationary 
epoch would explain why Qq is not very far from 1. 

For a fixed value of Qq, increasing to has the effect of shifting the curves up, as well as increasing somewhat the 
redshift at which the transition occurs. 

The two dashed curves in Fig. (|^) are shown for comparison, and represent c) dj, — c?L2(0.2, z) (horizontal dashed 
line), i.e. an open matter-dominated universe with fi TO = 0.2, and d) dz, = di,2(l,z) (lower dashed curve), i.e. a 
spatially flat matter-dominated universe. 



IX. THE AGE OF THE UNIVERSE 

As stated earlier, the only fundamental scale that enters into the effective action of the model presented here is to. 
Nevertheless, as we show now, the fit of our model to supernovae data predicts reasonable values for the age of the 
universe to. 

The relations derived in Section VII, leading up to Eqs. (|9^ ) and (|9^), give tj and to hi terms of to, ^o and Ho. 
For Ho = 65 km/ (s Mpc), m = 3.7 x 10~ 33 eV, and f^o = 0.4 (upper solid curve in Fig. |3|), we obtain 

tj = 5.66 x 10 9 years (113) 

t = 1.34 x 10 10 years. (114) 

For H = 65km/(sMpc), to = 3.2 x 10~ 33 eV and fl = 0.3, we obtain 

tj = 6.03 x 10 9 years (115) 

t = 1.33 x 10 10 years. (116) 

Therefore, in both cases a reasonable value of roughly 13 billion years is obtained for the age of the universe. More 
data at higher redshifts may lower the value of Qo m our model, which would further increase the age of the universe. 

Fig. |^ contains plots of the scale factor versus cosmic time for the two solid curves of Fig. [| In each case, the open 
matter-dominated universe that transits to the de Sitter phase is shown continued as a dashed curve, for comparison. 
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X. CONCLUSIONS 



In summary, we showed that a model in which a transition occurs from a matter-dominated to a de Sitter expansion, 
fits the SNe-Ia data. We call such a model, a transitionary universe. Also, we have proposed a free quantum field 
theory effective action, Eq. (|38|), in which such a transition evidently occurs. The existence of a particle of very 
low mass would cause the universe to make a transition from the matter-dominated to a new de Sitter stage. In our 
model, one can say that we are now observing the mass scale of this particle through the SNe-Ia data. 

Models involving interacting fields may also give a transitionary universe, and such models would be natural 
extensions of the free field model presented here as the simplest case. 

We emphasize once again that the solutions to the gravitational field equations we obtain, in particular the de 
Sitter solutions for £ < 0, exist without the necessity of a non-zero cosmological constant term in the effective action. 
Furthermore, these solutions are fairly insensitive to the presence of a cosmological constant term. In this manner, 
our model does not suffer from the problem of fine-tuning of the cosmological constant, which exists in mixed matter 
and vacuum energy models. 

The present matter density and the predicted age of the universe agree well with the current estimates. As in other 
models, a value of f^o not far from I may result from a period of early inflation. Further constraints on f2o and fn in 
our model could result from comparison to cosmic microwave background data, as well as from the time-temperature 
relationship during nucleosynthesis. We hope to carry out such a comparison in the future. 

Finally, we would like to mention that the -R-summed form of the effective action could have consequences for early 
universe cosmology as well. In particular, the existence of an imaginary term in the effective action, implying particle 
creation effects, could play a role in the exit from an inflationary universe. In future work, we plan to pursue these 
ideas as well as to carry out a dynamical calculation giving the details of the transition between the matter-dominated 
stage and the later de Sitter stage of the expansion. 
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APPENDIX A: DIVERGENCES IN THE EFFECTIVE ACTION AS M -> 0+ 

Consider the one-loop effective action in Eq. ( |30| ) below threshold, i.e. M 2 > 0. This action is then a real series 
which is divergent as M 2 — * 0. The part of the first term in this series, involving M 4 , vanishes in this limit. However, 
subsequent terms, involving / ; for I > 2, are all divergent as M 2 — * 0. All terms containing / 3 and higher orders 
correspond to an asymptotic series in inverse powers of M 2 and therefore do not give a valid expansion for small 
values of M 2 . The behavior of this expansion as M 2 — > is thus unphysical and does not seem to be a cause for 
concern^. However, the term f 2 ln(M 2 ) is also divergent in this limit and is physically required for renormalization 
of ultraviolet (UV) divergences and to obtain the trace anomaly. It is therefore of interest to examine in detail the 
divergent behavior of this term as M 2 — > 0. We will show here that this divergence, although infrared in nature, can 
be absorbed into the ultraviolet divergences of the theory by a procedure similar to the way infrared divergences at 
m = are handled in the usual one-loop effective action |l^] . To this end, we will consider the one- loop effective action 
(|30|), truncated up to the terms involving / 2 . We will work within the dimensional regularization scheme, which is 
more useful than the zeta function scheme in this context because the divergent terms are explicitly displayed. 

First, Eq. (Q) implies the following proper time representation of the effective action Jl0| , p^[ : 

1 r°° 

W l = —Tr dss^e-™ 11 

2 Jo 

= 1(^)2-^/2 y d D x J— g J°° ids ( ls )-l( 47ris )-^/2 e — (M 2 -»)p (XjXj j s ). (Al) 

The truncated one-loop effective action before renormalization is obtained, in four dimensions, by keeping the first 



three terms in a power series expansion of F in Eq. (15). These terms include UV-divergent contributions to the 
effective action arising from the behavior of the integrand near s = 0. All higher order terms are UV-finite. 
Thus we obtain for the truncated one-loop effective action 

W t \ un = ~(47r)-^V) 2 - B/2 / d °^[ ids(is)- D ^e-^ M2 -^(l + (^) 2 7 2 ), (A2) 

where we have set / 1 = without loss of generality. The above equation has a finite piece corresponding to the sum 
over all powers of s higher than 2 in the expansion of e~ ls t R in the integrand. Therefore Wtrun as defined above differs 
from what is usually regarded as the divergent part of the effective action by this finite piece. Here, we need to keep 
this extra piece in order to properly take the limit M 2 — > 0. 

In performing dimensional regularization about D = 4, it is convenient to define 26 = D — 4. One can perform the 



proper time integral in Eq. (A.2), to get 




<. un = (327T 5 *)- 1 / dW=S< A* 4 ( ~tt ) r(-2 -S)+f 2 (^P) r(-tf)| . (A3) 
Expanding the exponents and the Gamma functions about 5 = {D = 4), we get 

W t \ ua = (327T 2 )- 1 ! -2 + | -\HM 2 I^ 2 ) 

+ 12 (l~D 7 " ln ( M2 /^) }+°( D - 4 )- ( A4 ) 

The expression above already indicates that the divergence as M 2 — > may be absorbed into the UV-divergence as 
D —>■ 4 in the coefficient of / 2 . To see this explicitly, we may set M 2 = in Eq. (A.2) at the outset, and replace the 
upper limit of the s-integration by some large number T to regularize the integral (infrared regularization). We can 
then examine the divergent behavior as T — * oo. We therefore have 



8 It is possible that all higher order terms sum to give a finite contribution as M 2 — ► 0. An example of such a situation is the 
series (1 — s -1 ) -1 = J^^Lo x ~* wn ich is a valid expansion for x > 1. As x — + 0, every term on the right hand side is divergent, 
although the left hand side vanishes. 
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- d » xV =g j ids(is)-»r^ e -sz {1 + {is y f2) , 

which may be evaluated in terms of the incomplete Gamma functions 7(0?, x), 

WU = (327T 2 )- 1 ! d± X V- 9 {*- D/2 (e/n 2 ) D/2 l(-D/2,eT) 
+ r D ^ +2 J 2 (e/^) D / 2 - 2 7(2 - D/2, eT)} . 
We may now take the limit e — * of the above expression by using the power series expansion p5[ 

(-l) n x a+n 



This yields, after some simplification 



(32tt 2 ) 



2^-l 



*■ — ' n\(a + n) 

n— v ' 



f 2 



4-D 



+ ln( M 2 T) +0(D-A) 



(A5) 



(A6) 



(A7) 



(A8) 



As T — » 00, the first term in Eq. (A8) drops out, leaving behind a term proportional to / 2 which diverges logarith- 
mically in this limit. However, this divergence can be absorbed into the UV-divergence as D — > 4. 

We therefore find that the logarithmic divergence in the effective action in the "infrared" limit T — > 00 may be 
canceled by counterterms of the same geometric form as the ones introduced into the bare gravitational action to 
cancel the UV divergences as D — > 4, i.e. a counterterm proportional to f 2 is required here. The infrared problem in 
the partially summed form of the effective action is thus handled in a way similar to the infrared divergence in the 
usual effective action at m = 0, where a counterterm proportional to $1 is required. 

The analysis above may be carried out in a similar man ner even when we do not set M 2 — at the very beginning 
but rather let M 2 tend to zero from positive values. Eq. (A8) is recovered at the end of the calculation. 



APPENDIX B: VARIATIONS OF CURVATURE INVARIANTS 

Here we will list the variations of the curvature invariants that occur in the effective action in Eq. (^8|) . They are 
as follows: 



S( J d A x^R) = - J rfxy/^gSg^GP" 
J d i x^R 2 ) = J d A x^5 9llu l^g^R 2 - 2RR" V + 2R V » - 2. 9 ' ll/ Di?J 



- M- A IUR) + ^(2M~ 6 R^R^ - M'^R'^) + ^ w ln | A/ 2 /m 2 || 

5{ J d i x v ^g~R\n I A/ 2 /m 2 |) = J d i x y ^Sg^ {-G^ln | M 2 /m 2 \ -^M^RR"" 



M 4 



M 6 



S{ J d 4 x^R 2 \n I M 2 /m 2 |) = J d 4 x^g~5g^ jln | M 2 /m 2 | Qs"* ' R 2 - 2RR" U 

+ 2R ; ^) + £Ar 2 (-R 2 R" U + 6R ; "R U + 4RR' Uf *) - f.W l R{6R'"R' v 
+ RK' vfl ) + 2f'M- & R 2 R^R' v } 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 
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S( J d 4 x^f 2 ) = ± J d 4 x^Sg^ jsn^ - RW - \fTUR 
+2R^ VT R aT + l - g ^R aT R°- + 2R^ Tp R vaT P - l -g^R aTp xR° Tpx - AR^ff"^ 

S(J d 4 Xy/^g~nR In | M 2 /m 2 |) = J d 4 Xy/=g~ Sg^M' 2 {-2g^aUR + 2{UR)^ 

-^g^R' a R-,a + R'^R-v + lM~ 2 (v"(ni?) ;Q i? ;Q + 2</" y (ni?) 2 - AR^(oRy^ - 2drr^ 

+R^R a R, a + 2R a (g^D - VV)fi ;a + 2g^R. aP R af3 - 2R. a v R a ») 
+2l 2 M~ i (-R. a R' a (g^ v n - V»V V )R + AR^R a ^R a - 1//'" // " /.' : ] R. a ,) 

+ 6f M- e {g^R, a R a - R^R")} 

S( J d A x*/=g~R aP R afi In | M 2 /m 2 |) = J d i x^gSg^ jln | M 2 /m 2 \ [^g^ R afi R a[j - 3i^ a i^ 

+ uBT + ^ V UR - R» v + R av ^R a ^j + lAT 2 (-R^R af3 R afj + WUR - 2R a ^R, a ^ 

+g^R at3 R, afj + 2R. a R^ a + g^R- a R a - R^R U - 2R. a R a ^ - 2g^R aP UR a[i 
-2g^R a p, K R a ^, K + 2R aP '^R af) + 2R a ^R a ^) +f M~ 4 {-R^R, a R a - g^R a ?R ;a R.p 
+ 2R^R u)a R, cl + Ag^R al} R K R aP , K + g^R af} R af) UR - AR a p'^ R^ R afi - R af! R afj R^ 
+ 2fM- 6 R af3 R al3 (RVR» - // -7, , : ,,// : "i| 

S( J d A x^—gR a ^ s R aPlS In | M 2 /m 2 |) = J d 4 x^Sg^ jln | M 2 /m 2 \ {^g^ R a ^ 5 R a0lS 

-2R^ 5 R\ jS - AR^- f3a + lM- 2 (-R^R a ^ 5 R a0jS ~ 2g^ R a ^ b R a01& , R K 
-2g^R af3 ^R aPlS , K + 2R a ^R a ^ 5 + 2R a ^R a ^ - AR^ , a R, - AR^R, 0a ) 
+f M" 4 (Ag^R a ^ s R a ^R. K + g^ R a ^ s R aPl& UR - AR a ^ s R af3lS ^ R^ - R^ 8 R af3lS R^ 

+ AR Pl,a ^R, p R, a ) + 2£ 3 Af- 6 {-g^R, a R a + R»R V )} 
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FIG. 1. A plot of the LHS (bold-faced curve) and RHS (dashed line) of Eq. (49), as functions of y, for £ = 0.033, r = 10 
and A = 0. The slope of the dashed line increases as r decreases. 
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FIG. 2. A plot of the LHS (bold-faced curve) and RHS (dashed line) of Eq. (49), as functions of y, with £ = -0.03, r = 16 
and A = 0. As r decreases, the slope of the dashed line increases, and the intersection point is shifted closer to the value 
y = — £ . Recall that r — m 2 /m 2 P i and y = R/m 2 . 
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FIG. 3. A plot of the difference between apparent and absolute magnitudes, as functions of redshift z, normalized to an open 
universe with Oo = 0.2 and zero cosmological constant. The points with vertical error bars represent SNe-Ia data obtained 
from Ref.[5]. The two solid curves represent the values a) m = 3.7 x 1CP 33 eV and Qq = 0.4 (upper solid curve), and b) 
m = 3.2 x 10~ 33 eV and tto = 0.3 (lower solid curve). The horizontal dashed line represents an open universe with f2o = 0.2, 
and the dashed line curving downward represents a matter-dominated flat universe. Smaller values of f2o also would fit the 
data (see text after Eq. (112)). 
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FIG. 4. Two plots of the scale factor versus time for a spatially open model universe in which an initially spatially open 
matter-dominated cosmology evolves to a de Sitter solution. The parameters for the top model are m = 3.7 x 10~ 33 eV and 
Q.0 = 0.4, and for the bottom model, fn — 3.2 x 10~ 33 eV and flo = 0.3. The dashed curves represent a continuation of the 
open matter-dominated phase. 
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